Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 2 October 2009 



(MN WTe^ style file v2.2) 



Broad iron L-line and X-ray reverberation in 1H0707-495 



A. Zoghbi^*, A. C. Fabian^, P. Uttley^, G. Miniutti^, L. C. Gallo^, C. S. Reynolds^ 
J. M. Miller^ and G. Ponti^ 

^Institute of Astronomy, Madingley Road, Cambridge CBS OH A 

'^School of Physics and Astronomy, University of Southampton, Highfield, Southampton SO 17 IB J 
, ^LAEX, Centro de Astrobiologia (CSIC-INTA); LAEFF, PC Box 78, E-28691 Villanueva de la Canada, Madrid, Spain 
i ^Department of Astronomy and Physics, Saint Mary's University, Halifax, NS B3H 3C3, Canada 

^Department of Astronomy and the Center for Theory and Computation, University of Maryland, College Park, Maryland 207^2, USA 

^Department of Astronomy, University of Michigan, Ann Arbor, Michigan 48109, USA 
, APC Universite Paris 7 Denis Diderot, 75205 Paris Cedex 13, France 



2 October 2009 



ABSTRACT 

A detailed analysis of a long XMM-Newton observation of the Narrow Line Seyfert 
1 galaxy 1H0707-495 is presented, including spectral fitting, spectral variability and 
timing studies. The two main features in the spectrum are the drop at ~ 7 keVand 
a complex excess below 1 keV . These are well described by two broad, K and L, 
iron lines. Alternative models based on absorption, although they may fit the high 
energy drop, cannot account for the 1 keV complexity and the spectrum as a whole. 
Spectral variability shows that the spectrum is composed of at least two components, 
which are interpreted as a power-law dominating between 1-4 keV , and a reflection 
component outside this range. The high count rate at the iron L energies has enabled 
us to measure a significant soft lag of ~ 30 s between 0.3-1 and 1-4 keV , meaning 
that the direct hard emission leads the reflected emissions. We interpret the lag as a 
reverberation signal originating within a few gravitational radii of the black hole. 
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1 INTRODUCTION 



In a recent paper, Fabian et al. ('2009") reported on the dis- 
covery of the first broad iron L line in the Narrow Line 
Seyfert 1 galaxy 1H0707-495. The interpretation is that the 
line is broadened by relativistic effects very close to the black 
hole, where the line is thought to originate. 
This, along with the already established broad iron K 
line seen in many AGN, is a further confirmation that 
the observed X-ray spectra from suppermassive black holes 
in radiatively efficient AGN are emitt ed from with i n few 
gravitational radii of the black ho le (|Tanaka et al ] Il995l : 
iNandra et al. I I2OO7I . see lMilledlioOTi for a review). 

The standard picture is that matter is accreted in 
the form of an optically thick, geometrically thin disc 
(jShakura Svunvaevlll97 3') that radiates mainly in the UV. 
The usually observed power-law spectrum in X-rays, which 
is likely to be due to a Comptonising corona and/or the base 
of jet, is expected to illuminate the disc, producing a char- 
acteristic reflection spectrum that is observed alongside the 
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primary power-law (George & Fabian"'l99l"). The main fea- 
ture in the reflection spectrum is a fluorescent Fe Ka line at 
'^6.4 keV , which is expected to be broadened and smeared 
out if it is emitted close to the black hole. Other features 
are also expected to be seen in certain cir cumstances such 
as high iron abundance (^ Fabian et al.ll2009l , this work) , while 
a soft excess is naturally produced for a wide range of ioni- 
sation states. 

Alth ough alternative inter pretations of AGN spectra 
exist fe.g. lTurner Miller|[2009l for a review) , they generally 
fail to explain some variability properties seen in the data 
(see the discussion section). Spectral variability and timing 
analysis provide a completely independent way of studying 
the physics of emission. That, combined with the spectral 
fitting, can be used to distinguish between the different mod- 
els. In particular for the reflection interpretation, the model 
predicts a time delay, that in principle, is observable be- 
tween the prima ry emission and the reflected emission (e.g. 
IPoutanenI I2OOI). This has been seen for the first time in 
1H0707-495 (|Fabian et alll2009l ) and is discussed further in 
this paper. 

1H0707-495 {z = 0.0411, hereafter 1H0707), is a highly 
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Figure 1. Unfolded spectrum of 1H0707 in the 2008 XMM obser- 
vations (04) compared to the first observation (Ol). The spectra 
have been unfolded (using eeuf spec) against a power-law of index 
0. This removes the effect of the detector effective area and shows 
the main features in the spectrum. 04 represents the spectrum in 
the 'bright' state, and the Ol represents a typical spectrum for 
the 'faint' state. 



vari able Narrow Line Seyfert 1 (NLSl) galaxy (|Leighlvl 
Il99^ ). It has a steep spectrum and shows a strong soft ex- 
cess, and although these are c ommon among other NLSl, 
they are extreme in 1H0707 ('Vaugha n etaP Il999). The 
soft excess is traditionally interpreted as disc blackbody 
emission. However, when large samples of AGN are consid- 
ered, the measured temperature is higher than the temper- 
ature of a standard disc (Shakura & Syunyaev 1973), and 
independent of blackhole mass ( Gierlinski Pond 120041 ). 
It appears to be well explained by t he blend of relativis- 
tic ally smeared reflection continuum (|Crummv et al ] l2006l : 
see ISchurch Donell2008l for problems with alternative ab- 
sorption models). 

1H0707 gained interest after the discovery of a sharp 
drop in the spectrum at 7 keV , during the first XM M- 
Newton observation [01 hereafter] (|Boller et al.ll2002l ). A 
second XMM observation (02) sho wed a change in t he en- 
ergy of the edge from 7.1 to 7.5 keV (jCallo et al.ll2004 ). The 
drop was attributed to either the source bei ng partially cov- 
ered (|Boller et al.ll2002l : [Tanaka et al.ll2004l ). or it being the 



blue wing of a broad iron Ka line ( Fabian et al ] 120041 1. In 
2007, four short (40 ks each) observations were made with 
XMM (03 hereafter). A short report on a 4th much longer 
observation (04 he reafter) done in 2008 was presented in 
iFabian et al.l (|2009l l. In this paper, we present detailed anal- 
ysis of the results reported there. We also include analysis of 
the 03 data in some parts of this study which is organised 
as follows: observation and data reduction are presented in 
Sec. [21 spectral fitting and variability are presented in detail 
in Sec. [3] and |4] respectively. In Sec. O the lag calculation 
are shown, with all the results and implications discussed in 
SecM 



2 OBSERVATIONS & DATA REDUCTION 

1H0707-495 was observed for four consecutive orbits using 
XMM-Newton starting on 2008, January 29th, with a total 
exposure time of ^ 500 ks. The observations were made in 
the large window imaging mode. The data were analysed us- 
ing XMM Science Analysis System (SAS vS.O.O). The light 
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Figure 2. A plot of residuals in terms of sigma for the simple 
power-law + edge model fitted to the high energy band. 



curve was filtered to remove any background flares. No sig- 
nificant pileup was found in the patterns, which were se- 
lected using the condition pattern ^ 4 for the spectral 
analysis, and then relaxed when doing the timing. 

Source spectra were extracted from circular regions of 
radius of 35 arcsec centred on the source, and background 
spectra from regions on the same chip. The spectra were 
then grouped to bins with a minimum of 20 counts each. 
The response matrices were generated using rmfgen and 
ARFGEN in SAS. For the combined spectrum of the four or- 
bits, the events files were merged together before extracting 
the spectrum. 

The spectra show the presence of two instrumental lines 
between 8-9 keV . These appear as two clear emission lines 
in the background spectra, and as absorption lines in the 
background-subtracted source spectra. Their energies are 
8.0 and 8.9 keV respectively and can be identified as Cu- 
K lines originating from the electronics circuit board under 
the CCD (Katayama et al. 2004). We tried different back- 
ground regions to avoid the circuit board, but we found that 
for the best signal to noise ratio, it is optimal to get a large 
background region and include two absorption gaussians at 
8.05 and 8.91 keVin all fits. All results reported have been 
checked with spectra with cleaner background, but higher 
noise. 

Spectral f itting was performed using xspec vl2.5.0 
(jArnaudl [l996l ). All quoted errors on the model parame- 
ters correspond to a 90 per cent confidence level for one 
interesting parameter (i.e. a Ax^ = 2.71), and energies are 
given in the rest frame of the source (unless stated other- 
wise). The quoted abunda nces refer to solar abundances in 
lAnders k Grevessd (|l989l ). 



3 SPECTRAL ANALYSIS 

Fig.[T]shows the spectra from the 04 observation along with 
an earlier observation (01) for comparison. 04 shows a typ- 
ical spectrum when the source is in a 'bright' state, which is 
similar to 02, and three of the four 03 (not shown). The 01 
spectrum is typical when the source is faint which is similar 
to the first segment of 03. In what follows, the spectrum 
is explored in detail looking initially at hard energies (> 3 
keV), before including softer energies and looking at the 
whole spectrum. 



3.1 Hard Spectrum 

The main feature in the hard spectrum of 1H0707 is the 
sharp drop around 7 keV . To characterise the feature we 
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Figure 3. The spectrum of 1H0707 in the hard band (data 
points). The sohd hne is the model (a neutral edge from an out- 
flow with V = 0.03c) fltted to the data plus an emission line 
expected from simple physics given the edge depth. If the edge 
is due to ionised material, the fluoresence yield increases and the 
line should be stronger. 



fitted the 3-10 keV spectrum with a power-law and an edge. 
The fit was resonable (x^ = 313 for 258 degrees of freedom 
[d.o.f]). The edge energy is ^ = 7.3110.^6 keVwith a depth 
of T = 0.89t°:}3. The energy of the drop changed signif- 
icantly between earlier XMM observations, from 7.1 =b 0.1 
keVto 7.5 =b 0.1 keVfor the 01 and 02 observations respec- 
tively (jGallo et al.ll2004l ). The energy in 04 is between the 
two, and is more consistent with the 02 observation. The 
depth is also consistent within errors with that reported for 
02 (r = 0.841°:^^). 

Although the general shape of the drop is matched, few 
residuals clearly remain (see Fig. [2]), in particular the posi- 
tive excess at ~ 6.5 keVand the apparent absorption at '-^ 7 
keV . The former is possibly due to iron Ka emission. The 
latter might be caused by H-like iron absorption. The fit im- 
proves slightly if a smoothed edge is used instead, (Ax^ = 11 
for 1 less d.o.f), giving a width of 171^7° eV, consistent with 
01 and 02. If the spectral drop is interpreted as an iron K- 
edge, a corresponding emission line is expected. Adding a 
narrow gaussian line improves the fit slightly, it has energy 
of 6.731q"5q keV . The fit however, still leaves residuals be- 
low ~ 6.7 keV , which indicates that the line is broader. If a 
is allowed to vary, it has a value of a = 295lge^ i^n- 
isation states where resonant absorption is not important, a 
photon absorbed in the K-edge is followed by the emission 
of a fluorescent Une with a probability corresponding to the 
fluorescent yield. Its value is 0.34 for neutral iron and an ab- 
sorber that covers 47r solid angle as seen by the emitter. The 
observed line energy of 6.73 keVis inconsistent with it orig- 
inating in the same absorbing material (Eedge = 7.3 keV), 
regardless of ionisation or blue/redshift considerations. The 
90 per cent upper limit of the flux from the emission line 
in 1H0707 is 7.66 x 10"^^ erg s"^cm~^ , while the absorbed 
flux is 8.54 X 10""*^^ erg s~"'^cm~^ (the edge and line energies 
are now ^ 7.4 and 6.7 keV respectively) . This implies that 



I 

^ 0.01 



o 10"' 
5 

2.5 




-2.5 
-5 



\ \'iiV 

\ \f|L : 


Sk|. ' ' ' (b)' : 

^ 1 ILi - 
: 1 %^ 
1 % 

1 \ 1 — ! 1 1 't 






5 10 


- 1 - 

5 10 



Energy (keV) 

Figure 4. A plot of residuals in terms of sigma for the two models 
considered for the hard spectrum, (a): a partially covered power- 
law (green long-dashed) by an ionised absorber, and (b): a broad 
Ka line, (green long-dashed). For both plots, the dashed (red) 
line is power-law and the continuous (blue) line is the total. The 
data have been rebinned for display. 



< 9 per cent of the absorbed flux is re-emitted as fl ores- 
cence , while ^ 30 per cent is expected (Krolik & Kall manI 
Il987l l. Fig. [3] shows the line expected based on the absorbed 
flux if Att covering fraction is assumed. Although the ratio of 
expected to observed line flux might be reduced if the cov- 
ering fraction is less 20 per cent as implied by the data), 
it should be noted that this is just a lower limit given the 
ionisation assumed for the absorber. 

If iron abundance is included in the flt as a free param- 
eter (using zvfeabs in XSPEC), a lower limit (90 per cent 
confldence) on the abundance of ^ 18 x solar is found, while 
the best fit value is h igher, similar to earlier observations 
(e.g. boiler et al.ll2002l). Th e best fitting power- law index is 
r = 3.2 zb 0.1. Tanaka et al ] (|2004l ) used a broken power- law 
instead. That reduces the abundance slightly, but the fit in- 
dicates that the steepness introduced by the broken power- 
law is mainly driven by the deep drop, and not by an actual 
break in the power- law below 7 keV . A broken power- law 
fitted to energies below the drop does not fit better in a 
statistical sense than a simple power-law. 

To gain more physical insight into the K-shell absor ber, 
we u sed the XSTAR photoionisation code (K allman et al.l 
1 19961 ) to generate a grid of models to fit the data. We fixed 
the abundances at solar values except for iron which was left 
free, the illuminating fiux was a steep power-law of index 
3. Other fitting parameters include column density (A^h), 
ionization parameter (^ = AnF/n where F is the fiux at 
an absorber of density n) and covering fraction. If the iron 
abundance is fixed at solar, the fit is not good (x^ = 1.7 for 
258 d.o.f) and improves if it is allowed to vary, it has a lower 
limit of 14 X solar (panel (a) in Fig.|4|). At solar values, the 
high column density required to fit the drop gives too much 
continuum curvature below the edge. The best fit absorption 
model (which does not include any re-emitted components) 
has A^H 6 X lO^^cm"^, log^ ^ 3.0 and a partial covering 
fraction of ~ 0.7. This means that the absorber covers 20 
per cent of the sky as seen by the source, 70 per cent of it 
is in our line of sight. Given the high ionisation parameter, 
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Figure 5. A plot of residuals in terms of sigma for a model con- 
sisting of a partially covered power-law by a photoionised gas + 
disc blackbody (BASE), plus different phenomenological models, 
EDGE: Absorption edge, GAB: absorption gaussian, GA: emis- 
sion gaussian. 



the model also predicts a K/3 UTA fe.g. lPalmeri et alfcooj ). 
which is not seen in the data. If the absorber is allowed to be 
moving, an outflow velocity of ~ 0.01c is found (Ax^ = 12 
for 1 extra d.o.f). 

Alternatively, the spectral drop can be attributed to 
emission rather than absorption. The prominent emission 
feature in these energies would be that of iron Ka. The line 
has to be very broad to explain the drop and the excess emis- 
sion below ^ 6keV (a ^ 0.97 keV , E ^ 5.9keVand equiv- 
alent width of 1.2 keVif fitted with a simple g aussian). I f a 
relativistically broadened line is used (laor of lLaor|[l99ll in 
XSPEC), a good fit is found {xl = 1-03 for 255 d.o.f, see Fig. 
(31). The energy, inclination and inner radius of the line were 
found to be ^ = 6.19l° °I keV (rest frame), i = 53^1 de- 
grees and Tin = 1.56lQ Q9rg respectively, where rg = GM/c^ 
is the gravitational radius of the black hole. An equally good 
fit is found if the energy is fixed at 6.4 keVwith very small 
changes to the other parameters. The power-law index was 
r = 2.9 ± 0.1. No other significant emission or absorption 
feature is present in the spectrum (Fig. |4|. 

To see what effect the continuum would have on the 
broad line fit parameters, and to account for other emis- 
sion expected wi th iron line, we used t he self-consistent re- 
flection model of IRoss FabianI (|2005h reflionx. This was 
convolved with the kdblur kernel (in XSPEC) to account for 
relativist ic effects. The model has as free parameters the 
iron abundance, ionization ^ and the illuminating power- 
law index F (which is linked to the main power-law index 
in the fit), kdblur parameters include rin, i and q (where 
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Figure 6. The photoionisation absorption model fitted to the 
data. Main plot: the data and the model that under-predicts the 
emission at ~ 0.8 keV . Inset: a detailed plot of the model showing 
the UTA features due to iron. 



emissivity oc r ^), in addition to rout which was fixed at 400 
rg. The fit was very similar to that of the laor line. It has 
= 1.4 zb O.lrg, q = 8.9+^"^ and i = 45 ± 1 degrees. The 
iron abundance was found to have a low er limit (90 per cent ) 
of 7 X solar, similar to previous work (|Fabian et al.|[20o3 ). 
This indicates that most of the emission is coming from a 
small region very close the black hole. 



3.2 Soft Spectrum 

The spectrum of 1H0707 does not show the strong absorp- 
tion below ^ SkeVwhich would be expected given the ab- 
sorbing column density inferred from the deep Fe Ka edge. 
This means that, in the absorption interpretation, the ab- 
sorber only partially covers the source. The best fitting 
model consisting of a partially covered power-law was ex- 
tended down to 0.3 keV . A strong excess at soft energies 
(< 1 keV ) is present. T his feature is common in other NLSl 
fe.g lGierlihski Donei r2004). The standard picture is to at- 
tribute this to disc emission. However, the observed disc 
temperature (~ 150 keV) is higher than expected from a 
standard thin disc, and is independent of the bl ack hole mass 
(jGierlihski Donel[20o3 : ICrummv et al.ll2006h . Including a 
multi-colour disc blackbody improves the fit slightly (with 
temperature T = 0.14 keV ), however it leaves strong residu- 
als at ~ 1 keV (Fig.[5]top, this model is referred to hereafter 
as BASE) . The same residuals were found by iGallo et"al] 
(j2004 ) in this source, and by adding two gaussians (an 
emission and an absorption), the fit improves significantly 
(Ax^ - 7500 for 6 extra d.o.f, Fig.[5l). iiGallo et al.1 (|2004l ) in- 
terpret the absorption as being due to a blend of resonance 
lines, mostly from ionized L-shell iron, while the emission 
arises from an extended warm medium outside the line of 
sight. This seem rather arbitrary and, as suggested by the 
authors, a better explanation is required. 

The residuals are clearly not caused by absorption from 
OVII or OVIII, and including edges fixed at 0.74 and 0.87 
keV does not improve the fit at all. If the edge energy is al- 
lowed to vary it shifts to higher values but does not improve 
the fit by much (BASE+EDGE in Fig. [5)). The edge has en- 
ergy E = 1.011 ± 0.004 keVand r = 2.7 db 0.4. If this is a 
blue-shifted OVIII edge, the emitting material would have 
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Figure 7. A ratio plot of the data to a simple absorbed power- 
law. Galactic absorption was fixed at 5 x lO^'-'cm"^. 



to be moving at very high velocities 0.16c). A better fit 
is found if an absorption fine is used instead (Ax^ = 2181 
for the same number of d.o.f, with E = 1.163 =1= 0.004 and 
a = 0.0144 it 0.005). However, clear systematic residuals are 
left which indicates that the model is not consistent with 
the data (BASE+GAB in Fig. [5]). 

If the feature is modelled as an emission rather than 
absorption, a much better fit is found (x^ = 1495 for 1104 
d.o.f, BASE+GA in Fig. [5])- A fit with a gaussian emission 
line gives a line energy of ^ = 0.799 ±0.003 keV and a width 
of cr = 127 ± 6 eV. This is most likely due to iron L-shell 
lines. 

The super-solar abundances implied from the fits to the 
hard spectrum, can (or is expected to) have other possible 
effects on the spectrum. In particular, if the drop at ^ 7 
keV is caused by absorption, then one might expect to see a 
similar drop caused by the L-shell edge (^edge > 0.85 keV ). 
A simple phenomenological model consisting of a power-law, 
blackbody and two edges reproduces the general shape of the 
spectrum. This model is diflFerent from BASE+EDGE shown 
in Fig. [5]), where the full spectrum generated by XSTAR 
is used. This simple phenomenological model is unphysical, 
because it ignores the effect of the continuum, which needs 
to be modelled properly in any absorption scenario. 

The fact that an absorber producing the K-edge has to 
be partially covering the source imply that the two features 
at 7 and 1 keV cannot be due the same absorbing cloud. The 
best fitting (partially covered) XSTAR model of Sec. l3.1l was 
extended to lower energies, a second ionised absorber was 
added to try to account for the ~ 1 keV feature. A black- 
body component was also included to account for the soft 
excess. The parameters of the second absorber were allowed 
to vary. No good fit was found. The main reason is that the 
model tries to fit the apparent drop below 1 keV with an L- 
shell iron edge, requiring high abundances, however, as Fig. 
[6] shows, this predicts a drop between 0.7 and .9 keVdue to 
the M -shell unresolved transition array (UTA) (|Fabian et al.l 
I2OO9I I. This is a blend of numerous absorption lines arising 
from the photoexcitation of ions Fel - FeXVI m ainly pro- 
duced by 2p — 3d transitions (Behar et al.ll200lh . It should 
be noted that OVII and OVIII edges are also predicted by 
the model but not seen in the data, particularly for smaller 
values of the iron abundances. iNicastro et aU (| 19991 ) showed 
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Figure 8. The best fitting reflection model with data to model 
ratio. Long-dashed (red) line is the power-law, dot-dashed (green) 
is the blackbody and dotted (blue) line is the blurred reflection 
component, the black line is the total. The reflection fraction 
(ratio of the reflection to power-law fluxes) is 1.3. 

that similar features in an ASCA observation of IRAS 13224- 
3809 can be reproduced by a warm absorber that is illumi- 
nated by a steep power-law. This, however, predicts narrow 
absorption lines whic h are not present i n the RGS data of 
1H0707 (see Fig. 3 in lFabian et"al]l2009l l. 

The conclusion of this analysis is that the spec- 
trum around 1 ke V can not be described by absorption. 
Fabian et al. 1 (|2009") has pointed out that this spectrum can 
be due to a broad Fe-L line. This is expected if the high 
energy drop is interpreted as the blue wing of a broad Ka 
line. Their ratio plot (their Fig. 9, which is reproduced here 
in Fig. [T]) of the data to a simple model consisting of a 
power- law with galactic absorption show clearly two main 
broad line-like residuals. If, as a phenomenological model, 
two Laor fines are added to the fit, a good description is 
found (though not statistically acceptable) . When the fit pa- 
rameters of the two Laor lines are linked, the fit parameters 
are: Hn = 1.82 zb 0.04 rg and IncL = (5li;;^)13. The emis- 
sivity index is 5.3 =b 0.1. The energies of the lines are ^0.91 
and ^ 6.0 keV (rest-frame) for the L and K lines respec- 
tively. These lines are not exactly at the expected energies 
because the continuum is not modelled properly, and this 
emphasises the point that these lines cannot be fitted sepa- 
rately, and need to be modelled as part of the whole reflec- 
tion continuum. If the line parameters between the two lines 
are allowed to vary independently, they still give similar and 
consistent results (K-line: nn = l-4lo".2 ^g? ^^c/. = 57 =b 5°, 
L-line: nn = 1.41±0.05 rg, Incl. = 55.6itl.0°, the emissivity 
index is 5.8 di 0.2 in both cases). 

To account for the whole reflection spectrum expected 
to accompany the two broad lines, we again used the 
relativist ically-blurred self- con sistent reflection model RE- 
FLIONX (|Ross FabianI l2005l l . The best fitt ing model is 
shown in Fig. [8] (this is similar to Fig. 10 in iFabian et al.l 
|2009' shown there for the first orbit only) . The model gives a 
very good description of the data. A low temperature black 

^ The fit paramete rs are slightly different from those reported in 
I Fabian et al.l (|2009l ) because only the flrst orbit of the data was 
used in getting their parameters, whereas all four orbits are used 
in this work. 
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Figure 9. The light curve of the current observation compared to previous XMM observations. Ol, 02 and 03 refer to the first, second 
and third observations respectively. 04 is the current observation. Ol was made in full frame mode, while 02, 03 and 04 were made in 
large window mode. The vertical dashed lines separate the different segments of the observation. 



body was required to fit a small excess below ~ 0.5 keV . 
This is different from the blackbody generally used w hich 
have higher temperature {r^ 0.15 keV , iGierlinski Pond 
I2OO4I ) . The temperature of the disc blackbody in our fits is 
kT ~ 50 eV, which is more consistent with what is exp ected 
from a thin accretion disc fehakura Svu nvaevllTQTih . 

The fit parameters for the blurred reflection model were 
consistent with the two- /aor flt. The inner radius was Tin = 



0.8 



0.9 



1.1 



1.2 



1.23Z 



with a disc inclination of Incl. — 58.5 



+0.8 
0.7 



The emissivity index was 6.6_]^j^. The combination of small 
inner radius and emissivity index is an indication that most 
of the emission originates very close to the black hole, and 
that explains the very broad iron lines. The power-law has 
a steep spectrum of ~ 3.2, this is in line with the estab- 
lished fact NLSl h ave steeper than usual X-ray spectra 
terand t et al.|[l997l l . The residuals in the flt between 2 and 5 
keV can be explained by ionisation changes of the spectrum. 
The data fltted is the time average of four orbits, and trying 
to parametrise a highly variable spectrum with a single ion- 
ization model might not be physical. If an extra reflection 
component is added, with all parameters linked to the flrst 
except for the ionisation (and normalisation), the excess dis- 
appears. For the small residuals at ^ 1 keV , most of them 
are due to ionisation changes, in addition to elements other 
than iron having non- solar values. The process responsible 
for producing the high iron abundances, inferred from the 
K and L lines, is also expected to affect the abundances of 
other elements. 



4 SPECTRAL VARIABILITY 

The source has shown signiflcant variabi Uty in the past. In 
particular, the flrst observation of 2000 (|Boller et al.ll2003 l 
seems to have caught the source in a low state compared 
to the second observation (|Gallo et al ]|2004). The source in 
the current observation is in a similar flux state to that of 
02, which is clear from the flt parameters above, and the 
light curve shown in Fig (9] 
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Figure 10. Time-resolved spectra for lOks segments. Top: a 2D 
image of the ratio of the spectrum in each segment to the average 
spectrum, as a function of energy and time, a typical error in 
these points is ±0.03. Bottom: the corresponding count rate in 
each segment (0.3-10 keV). 



4.1 Time- resolved spectra 

The spectral fltting presented in Sec. [3]showed that a reflec- 
tion spectrum originating close to the black hole can flt the 
data very well, with alternative absorption models failing to 
account for the complexities at ^ 1 keV . In this section we 
investigate the variability of the spectrum with time. This 
can provide more constraints on the spectral models. 

The 04 data have been divided into segments of 10-ks 
length. A ratio of each spectrum to the average was calcu- 
lated. This will highlight the variability pattern in the spec- 
trum on the time-scale probed. This was achieved by fltting 
the time-averaged spectrum with a multi-spline model to 
get the best description of the data, and then fltting it to 
the segments allowing only a multiplicative constant to vary. 
The spectra in each segment were binned using grppha so 
that each bin has a minimum of 20 counts per bin (with 
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Figure 11. Time-resolved spectra for 20ks segments. The data 
is presented in a similar way to that in Fig. [To] The data image 
have been slightly smoothed to show the patterns more clearly. 



similar bins between segments), then each ten energy bins 
were grouped together so the errors are further reduced. 

The results are shown in Fig. [101 alongside the cor- 
responding count rate from each segment for E < 3.5 
keV where the noise effect of noise is minimum. The first 
point to notice in the variability is the apparent line di- 
viding the spectrum at ~ 1 keV . The variability patterns 
below and above 1 keV appear anti-correlated. This may be 
an indication that there is spectral pivoting at 1 keV . 

The other noticeable pattern is the fact that residu- 
als above 1 keV match the total count rate very well (Fig. 
nop . As the total flux gets higher, the residuals above 1 
keV increase and those below get smaller. The opposite is 
also true. 

To investigate the possibility of spectral pivoting at 1 
keV , we extended the analysis to include higher energy bins. 
This was not possible with 10 ks segments, so we used 20 ks 
instead. Fig. [TT] shows the results. The energy now extends 
to ~ 8 keV . The data were binned in a similar way to Fig. 
[To] except for grouping every 2 instead of 10 energy bins. 

The figure shows similar patterns to those in Fig.[Tni but 
now the variability pattern above 1 keVdoes not continue 
up to higher energies, and there appear to be another cut at 
~ 5 keV . Also, there seem to be indications that variability 
above ~ 5 keV is similar to that below 1 keV . 

The fact that the variability pattern does not extend 
above 5 keVis inconsistent with spectral pivoting. It can 
rather be explained if the spectrum is composed of at 
least two components, one dominating between 1 and 5 
keV (hereafter C<i-5>), which is highly correlated with the 
total flux of the source. The other component dominates 
outside the ^1 — 5 keV range (hereafter C>i-5<). 

This model-independent description of the variability is 
fully consistent with the reflection model presented in Fig. 
[8] A power-law component (PLC) dominates between 1 — 5 
keV and a reflect ion- dominated component (RDC) is present 
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Figure 12. Left: Flux-resolved spectra for 10 flux segments, each 
column in the image is a ratio of the segment spectrum to the 
average. The data image have been slightly smoothed to show the 
patterns more clearly. The error in points is of order 0.03. Right: 
The same as in the left panel but obtained from simulations in 
which the power-law varies up and down by 20 per cent keeping 
the reflection constant (see. Fig. [8]). 



outside this ra nge. This picture is similar to that reported for 
MCG-6-30-15 (jFabian Vaughanlliooi : I Vaughan FabianI 
E)04l l 



4.2 Flux-resolved spectra 

Another way to exp lore the variability is t hrough flux- 
resolved spectra fe.g. IVaughan FabianI l2004l ). We carried 
out similar analysis to that in Sec. 14.11 but now the spectra 
are extracted from flux segments instead of time segments. 
The light curve was divided into ten flux bins with each hav- 
ing equal number of points. The spectra from each segments 
was then compared to the average spectrum. 
Fig. [12] (left) shows the result. At low fluxes the spectrum 
shows an excess (compared to the average) below 1 keV and 
above 5 keV . As the flux increases the excess is shifted, 
and most of it is between 1 a nd 5 keV . This, as was 
pointed out in Sec. 14.11 (see also iFabian Vaughanl l2003l : 
IVaughan Fabianll 20041 ). can be easily explained if the spec- 
trum if composed of two components (C<i-5> and C>i-5<), 
and seem to be consistent with a power-law component 
(PLC) between 1 and 5 keV and a reflection-dominated out- 
side the range, as described in Sec. [3] 

The two-component interpretation of the spectral 
variability has be en su ggested by d i fferen t authors for 
Seyfert galaxies (iFabian Vaulchanl l2003l : iTavlor et al.l 
I2OO3I : IVaughan FabianI 1 2004l l. 1H0707 seems to show sim- 
ilar patterns. This simple picture is also consistent with the 
best fitting model shown in Fig. [8] For example. Fig. [12] 
(right) shows the result of applying the same analysis pro- 
cedure of Sec. l4.2l to a simulated set of spectra, produced by 
taking the best fitting model of Fig. [8] and generating differ- 
ent spectra by changing the normalisation of the power-law 
randomly by 20 per cent, and it shows clearly that the pat- 
terns are reproduced very well. 

Another representation of the spectral changes as a 
function of flux is shown in Fig. [13] It shows a ratio of the 
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Figure 13. Ratio of the spectra in each flux segment to an ab- 
sorbed power-law of index 3. Flux increases up. The y scale is 
shown for the flrst plot, others are plotted on the same scale. The 
horizontal line with each plot represents the ratio of 1. 



spectrum in each flux segment to an absorbed power-law of 
index 3. 

What is apparent is that the spectral drops, both at 1 
and 7 keV become weaker as the flux increases. This, com- 
bined with the earlier conclusion that the spectrum below 1 
keVand above 5 keV represents the same component, point 
to one conclusion: C<i-5> is well correlated with the total 
flux, while C>i-5< is more constant, the contrast between 
the two components causes the observed effect. 



4.3 Difference Spectrum the Nature of the 
Spectral Components 

We showed earlier, using flux and time-resolved spectra, that 
the spectrum appears to be composed of two components. 
To find their nature, we plot th e dif ference spectrum in Fig. 
fT4l ^similar to Fig. 7 in lFabian et al.ll2009l ). This is produced 
by taking the difference between two flux segments (above 
and below the mean count rate) . This way, if the spectrum is 
composed of two components, a constant and a variable, the 
constant part is subtracted and only the variable component 
is left (jFabian Vaughanll2QQ3l ). 

The variable component cannot be fitted with a single 
power- law (x^ = 2.9), however if the power-law is fitted only 
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Figure 14. The difference spectrum between high and low states, 
deflned as that above and below the mean respectively, along with 
and absorbed power-law flt and the residuals. This represent the 
variable component. 
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Figure 15. A plot of the reflection flux vs power-law flux (both in 
0.3-10 keV ) changes in the 10 ks segments. The best fltting model 
in Fig. Owas fltted to each segment, allowing the normalisations, 
photon index and ionisation to change. 



above 1.5 keV , a very good fit is found with an index of 3.2 
(X^ = 218 for 270 d.o.f, see Fig.[Tll). 

Although interpreting the physical origin of the variable 
component might not be straight forward, the constant com- 
ponent is consistent with the C>i-5< discussed earlier, and 
can be easily identified as a reflection dominated component 
(RDC), characterised by the two broad K and L lines and 
the accompanying reflection continuum. 

The variable component is well-fitted with a power- 
law above ^ 1.5 keV , an d is similar to that fou nd in 
oth er AGN (MCG-6- 3Q-15 Vaughan h FabianI l2004l . Mrk 
766 iMiller et al.1 l2007l . although both sources show signa- 
tures of warm absorption, unlike the case here). We interpret 
the variability below 1 keVas being due ionisation changes 
(iFabian et al.1 [2009) . The spectrum, as the best fit shows, 
is dominated at these energies by reflection, mainly iron 
L emission and emission from other elements, and as the 
power- law changes, we would expect the reflector's ionisa- 
tion to change too. 

To investigate this interpretation further, we fitted each 
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of the time- resolved spectra discussed in Sec. 14.11 individu- 
ally with the best fitting model shown in Fig. [8] The nor- 
malisations of the power-law and reflection components were 
allowed to vary, as well as the power-law index and the ion- 
isation of the reflection. Other parameters are not expected 
to vary on the time-scales probed here so they were fixed. 
To increase the ranges of parameters, we also included data 
from the two earlier XMM observations of the source (01 
and 02) . Fig. [15] shows a plot of the flux changes of the 
power- law and reflection components inferred from the fits 
to 10-ks segments. The figure shows that the two compo- 
nents are correlated, particularly if the 01 low state data 
(red triangles) are also included. It is also apparent that de- 
spite this correlation, most of the variability is driven by the 
power- law component. Two orders of magnitude change in 
the power-law corresponds to a factor of ^ 3 change in the 
reflection component. This is generally similar to the pattern 
seen in MCG-6-30-15 (Revnolds et al. 2004), where at low 
fluxes the reflection component is correlated with the power- 
law, while at higher fluxes the reflection component is more 
consistent with a constant. Also, we note that there appears 
to be a hysteresis behaviour in the relation at high fluxes 
(mainly in 04). If real (hard to say), it could be related to 
the geometry of emitting region. 

In a reflection scenario, the reflected radiation is ex- 
pected to respond to the primary power-law variations, and 
a positive correlation is expected. The fact that this is not 
the case at high fluxes is usually explained by light bendin g 
effects (IFabian k VaughanI [2OO3I : iMiniutti k FabianI l2004l ) . 
The variability is driven mainly by a primary source moving 
in the strong gravity field. This can produce up to an order of 
magnitude flux variations with small changes in the reflected 
flux. The decreasing strength of the featu res at 1 and 7 
keVis similar to that seen in MCG-6-30-15 (jVaughan et al.l 
I2OO3 II and is in line with the pre dictions of the light bending 
model (Mini utti Fabian|[20o3 ). 

The other interesting result is related to the ionisation 
changes. Fig.ll6lshows the variations in ionisation parameter 
^ along with the power-law flux, and shows that the ionisa- 
tion parameters changes as the power-law flux changes. In 
agreement with the interpretation of the difference spectrum 
of Fig. [m (see Sec.[6]for more interpreting this in the light of 
the light bending model). The ionisation changes of the re- 
flector between 10-100 (erg cm s~^) causes the shape of the 
reflection spectrum (that is possibly otherwise constant) at 
^ 1 keV to change, it is mainly caused by changes in the iron 
L emission, and also possible changes from other elements, 
and this is responsible for the 'bump' in Fig. [M] 



5 TIMING ANALYSIS 

To further understand the properties of the source, we focus 
in this section on the timing properties to establish whether 
or not it is consistent with our interpretation of the spectrum 
and spectral variability. 

One key prediction of the reflection interpretation of the 
spectrum is a time delay (lag) between the primary illumi- 
nating power- law and the emission from the reflector acting 
as a mirror. 
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Figure 16. Ionisation of the reflector vs power-law flux (0.3-10 
keV) in the 10 ks segments fits (see also Fig. llSland text). 



5.1 Coherence 

The coherence of two light curves is a measure of how cor- 
related they are, or how much of one light curve can be 
predicted from the other. Mathematically, it is defined as: 

\{C{f))f 



lif) 



(1) 



(I5(/)P>(IH(/)P) 

where S{f) and H(f) denotes the Fourier transforms of 
the two light curves under study. C{f) = {S*{f)H{f)) is 
the cross spectrum, and * denotes the complex conjugate. 
The angled brackets indicates the average over several light 
curves and/or frequency points. 

The coherence takes a value of 1 if the two light curves 
are perfectly coherent (for details on the meaning of coher- 
ence and how it is calculated, see Vaugha ^rA"N owak| Il997l: 
iNowak et al.lll999l : F^ughan et al.ll2003l ). 

Fig. [T71 shows the coherence function 7^(/) between the 
two energy bands 0.3-1.0 (soft) and 1.0-4.0 keV(hard). To 
calculate this, background- subtracted light curves were ex- 
tracted in the two bands, then divided into 4 segments of 
about 100 ks length each. The Fourier transform of each 
segment was then taken and the coher ence and its errors 
were calculated following equation 8 in IVaughan Nowa3 
((1993) . The frequency bins were constructed so that the bin 
size is equal to 1.4 times the frequency value. The figure 
shows that the two Ught curves are highly coherent 1) 
for almost all of the frequency range of interest. The coher- 
ence drops sharply at 5 x 10~^ Hz where the noise starts 
to dominate. 

To confirm this we simulated 1000 light curve pairs 
(each representing the two energy bands), that have the 
same statistical properties as the observed data (see Sec. 15.31 
on details about the simulations). We followed the method of 
Tim mer k Koenig) (| 19951 ) to generate red- noise light curves 
that have the same mean and variance as the data, then 
applied Poisson noise and the same windowing and binning 
functions as that used for the real data. The simulated light 
curves are perfectly coherent and the only source of devi- 
ation from 1 in the coherence function is due to Poisson 
noise. As can be seen is Fig. 1171 where the blue continuous 
line is the median of the coherence measured from 1000 light 
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Figure 17. Black dots: Coherence as a function of Fourier fre- 
quency (inverse of time-scale) between the two energy bands 0.3- 
1.0 (soft) and 1.0-4.0 keV(hard). The blue line is the mean of 
the coherence of 1000 simulated light curve pairs that have the 
same properties as the observed data. The dashed lines represent 
the 95 per cent confidence levels. The frequency bins were con- 
structed so that the bin size is equal to 1.4 times the frequency 
value 



curves, the data is fully consistent with the simulated light 
curves, which imply that the light curves of interest are high 
coherent up to ^ 5 x 10~^ Hz. 



5.2 Time Lags 

Fig. [T8l shows the time lag between the 0.3-1.0 and 1.0-4.0 
keV energy bands as a function of Fourier frequency (similar 
to Fig. 3 in iFabian et al. 1[2009). A similar procedure to the 
coh erence was fo l lowed . The lag t(/) was calculated follow- 
ing |N^^etd] (|l99i) with 



rif) 



m _ arg[C(f)] 



27r/ 



27r/ 



(2) 



where C{f) is again the cross spectrum, and arg[C(f)] is the 
argument of the complex number C{f). The sign convention 
here means that a positive lag indicates the soft flux changes 
before the hard flux. The erro rs in the la^ were c alculated 
again following equation 16 in lNowak et al.l (|l999 h. and are 
investigated using Monte Carlo simulations in Sec. 15.31 

The figure shows that below ~ 5 x 10~^ Hz, the soft flux 
leads the hard flux by about 150 seconds, the lag appears to 
turn over below ^ 10~^ Hz. At ~ 6 x 10~^ Hz, the lag turns 
negative, and variations in the hard band now lead. However, 
above ~ 5 x 10~^ Hz, the coherence in the light curves is 
lost fFig. ll7|) and the lag is dominated by Poisson noise. The 
shape of the lag function cannot be characterised by a single 
power- law as it is commonly done with other sources (e . 
IPapada kis et al.1 l200ll : IVaughan et~aD l2003l : lArevalo et al 
2006). However, if fitted only to the range [2x 10"'^-3x 10"^] 
Hz, a power- law gives a good description with the best fit 
function being r(/) = 0.34/-°-^^ - 87. 

The lag measurements could be affected by artefacts 
and biases in the Fourier calculations, and Poisson noise 
could also have a contribution. To investigate the signifi- 
cance of the lag measurements, the best procedure is to use 
Monte Carlo simulations, and that is the topic of the follow- 
ing section. 
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Figure 18. Time lag as a function of Fourier frequency (inverse 
of time-scale) between the two energy bands 0.3-1.0 (soft) and 
1.0-4.0 keV (hard). The frequency binning is similar to that of 
Fig. [17] A positive lag (mainly below ~ 5 x 10~^ Hz) indicates 
that the hard flux lags behind (i.e. changes after) the soft flux. 



5.3 Lag Significance 

The idea here is to generate light curves with the same sta- 
tistical properties as the data, and then impose a defined 
lag and see how well it can be recovered. 

The light curves were generated fro m red-noise power 
spectr a (PSD) following the method of iTimmer Koenid 
(| 19951 ). The PSD of the data was calculated using light 
curves from the whole energy band, and was fitted with a 
broken power-law. The best fitting model had an index of 
2.18 ±0.09 at high frequencies breaking to 0.8 ±0.5 at a fre- 
quency of /break = 1.41q 4 X 10~^ Hz (if the low frequency 
index is fixed at 1, the break is 1.61q q x 10^ Hz). This best 
fitting model was used as the underlying PSD in generat- 
ing the simulated light curves. Using energy-resolved PSDs 
instead of the total PSD does not change the results of the 
simulation. 

One thousand light curve pairs were generated in each 
experiment, where for each light curve pair, a random re- 
alisation of the Fourier transform for one light curve was 
generated. The second light curve Fourier transform was 
produced from this by adding a frequency-dependent phase 
to the Fourier transform. They are then inverse Fourier- 
transformed to yield two light curves that have the desired 
lag function between them. The light curves were scaled so 
they have the same mean and variance as that of the data 
before adding Poisson noise. To account for any bias in the 
data that usually are associated with Fourier transforms, 
mainly windowing and aliasing effects, the generated light 
curves were longer than the data, they were then resampled 
to match the data. 

Fig. [19] shows the result of averaging 1000 lags mea- 
sured for three sets of light curves, with input lags of 0, 
—30 and ±100 seconds, so they roughly cover the range of 
lags observed. The error bars represent the la spread in the 
measured lags from different light curves. The lags can be 
recovered to a very good degree in the middle frequencies 
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Figure 19. Time lag as a function of Fourier frequency for three 
sets of simulated light curves. The time lag was calculated from 
each light curve and then averaged. The error bars are the stan- 
dard deviation in the spread of the measured lags. Blue squares 
are for a lag of 100 s, red diamonds are for s and black circles 
are for -30 s laff. 
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5 X 10~^ — 5 X 10""^ Hz, with the probability of measuring 
a lag that is off by more than 15s from the true value is 
less than 5 per cent. 

At higher frequencies, as we saw in Fig. [m the light 
curves coherence drops below 1 sharply, and the noise dom- 
inates the data. For the 100 s lag, the sharp drop in the lag 
is caused by the fact that the lag is comparable to the fre- 
quency over which it is measured (i.e. the lag is of the same 
order as the time-scale it occurs). The drop happens when 
the phase suddenly jumps from tt to — tt, which for 100 s 
lag occurs at / = (j)/{2'KT) = 5 x 10~^ Hz. This is also the 
reason the lag tends to zero at high frequency. Phase flipping 
between tt and — tt makes the lag averages to zero. At lower 
frequencies, the spread in the measured lag (or equivalent ly, 
the errors in Fig. I18|) becomes large, and that is because 
the time-scale probed now are comparable to the segments 
size over which the Fourier transform was performed. Also, 
it should be noted that spread in the lag measurements is 
fully consistent with the error formula in Nowak et al... (.1999, ) 
(their equation 17). 

To see if the shape of the lag function itself (other than 
constant) has any effect on the measurements, we repeated 
the simulation for three lag functions: a power-law of the 
form r(/) = 0.34/-°-^^ - 87 which was found to fit the cen- 
tral frequencies of the measured lag, a step function that 
changes from +150 s to —30 s at a frequency of 4 x 10~^ Hz, 
and finally an interpolated smooth function describing the 
measured lag. The result are shown in Fig. [20l The lag func- 
tions are clearly well recovered, with the earlier mentioned 
trends visible at very high and low frequencies, with the lag 
at the very low end slightly under estimated. 

The main conclusion is that for most of the frequency 
range of interest, the spread in the simulated light curve lags 
caused by the noise is small and the lags measured in Fig. [18] 
are significant. In particular the soft lag which is the first to 



Figure 20. Time lag as a function of Fourier frequency for 
three sets of simulated light curves. The input time lag func- 
tions are (from top to bottom): a power-law of the form t(/) = 
0.34/~*-'-^^ — 87, a step function that changes from +150 s to —30 
s at a frequency of 4 x 10~^ Hz, and an interpolated function 
describing the measured lag of Fig. [18] The points are the mean 
of the measured lag from 1000 light curve pairs. The error bars 
represent the one sigma spread. The input lag function is plotted 
as a dotted line in each case. 

be se en in active galaxies ([Fabian et al ."2009'. see iGallo et al.l 
l2004l for a possible soft lag in IRAS 13224-3809). 



5.4 Alternative Approach: Time domain 

A significant lag in the data, in particular the soft lag, would 
in principle be visible in the time domain light curves. The 
problem however, is apparent from Fig. [181 and that is, a 
raw light curve would, if anything, show the average lag 
over all time-scales. In order to get any useful lags in the 
time domain, the time-scale (or a frequency range) has to be 
selected first. This is achieved by smoothing the light curve 
by two gaussians of slightly different sizes, then subtracting 
them, so that only variations on the time-scale of interest are 
left. This is the 1-dimensional equivalent of unsharp masking 
in image processing. The shape of the smoothing function 
is unimportant for the study presented here, and several 
functions were tested and they do not alter the conclusions 
drawr0. 

Fig. [2T] shows two segments of the light curves at the en- 
ergies of interest (0.3-1.0 and 1.0-4.0 keV ). In the top panel, 
the time-scale corresponds to the lowest negative point in 

^ Some functions produce biases in the lag that can be calibrated 
using Monte Carlo simulations 
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Figure 21. Two segments of the light curves at 0.3-1.0 (red dotted line) and 1.0-4.0 keV (blue continuous line). Top: variations on 
time-scale of 700s, corresponding to the lowest negative point in Fig. [18] Bottom: variations on time-scale of 5000s, corresponding to the 
highest positive point in Fig. 1181 Notice the time-scale on the x-axis. 
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Figure 22. The average CCF (see method (i) for details), for two 
time-scales, 700s and 5000s, corresponding to the minimum and 
maximum in Fig. 1181 The vertical dotted line is the zero lag line. 
This clearly shows that the CCF is distorted towards negative 
lags for the 700 s case, and towards positive lags for the 5000 s 
time-scale. Notice the difference in x-axis scaling. 

Fig. nil and it shows how the dotted line (0.3-1.0 keV) lags 
behind the continuous curve (1.0-4.0 keV), and this is the 
negative soft lag seen in Fig. [18] The trend is reversed in the 
bottom panel when a different time-scale is selected. Now 
the dotted line leads, and this corresponds to the positive 
lag at low temporal frequencies. 

To quantify this lag in a systematic way, we employed 
two methods that measure the lag as a function of time-scale. 
These are essentially equivalent to measuring it through 
phase lags using Fourier transforms. 

(i) Correlation: After selecting the time-scale of interest 
(e.g. Fig. [2T]) . the light curve is split into segments whose 
length is several times the time-scale (e.g. for a time-scale 
of 700 s, the length of the segments is 7000 s). The cross 
correlation function (CCF) between light curves in the two 
energy bands (0.3-1.0 and 1.0-4.0 keV) is then calculated 
for each segments. The CCFs from the segments are aver- 
aged, and the lag is measured through the maximum of the 
average CCF (see for example Fig. [22]) . 
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Figure 23. The same as Fig. 1181 with the result of the two 
time domain methods over-plotted. Black dots represent the lag 
measured using the phase difference in the Fourier domain. The 
shaded areas represent the la error for the lag measured using 
the methods discussed in Sec. 15.41 The lag was calculated for 100 
time-scales between 10^ and lO'^ s, with the errors coming from 
simulations of 50 light curve pair for each point. 

(ii) 'Flare' Stacking: In this method, again after selecting 
the time-scale in the light curves of the two bands (0.3-1.0 
and 1.0-4.0 keV) using a smooth function (e.g. Fig. [21]), all 
peaks in one of the two light curves (which is used as a refer- 
ence) are identified. For each peak, the corresponding '/?are' 
is identified by selecting the points either side of the peak. 
These segments of the reference light curve are matched by 
their time equivalent from the second light curve (not nec- 
essarily representing flares). The lag is found by comparing 
the average of all flares (.i.e an 'average flare'), with the cor- 
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responding (average) segment from the second light curve 
(this is to say, on average, what does a flare in one band 
correspond to in the other band). 

In order to obtain an estimate of the errors of the lag 
measured using the two previous methods, we again used 
Monte Carlo simulations, where for each measured lag, 50 
light curve pairs that have lag equivalent to the measure- 
ments are generated, the same procedure of Sec. l5.3l is used, 
the error is taken as the la spread in the distribution of the 
measured lags. 

Fig. 1231 shows the result of doing the calculations using the 
two methods for 100 time-scales between 10^ and 10^ sec- 
onds, the shaded areas represent the errors taken as the la 
spread in the lag calculation from simulated 50 light curve 
pairs for each point. This method was tested extensively on 
simulated data with known lags, and was found to recover 
the lags very well. This figure shows a very good match be- 
tween the methods and further emphasises the significance 
of the lag measurement. 



6 DISCUSSION 

A detailed analysis of a long XMM observation of 1H0707 
has been presented. It is clear that the source during the 
present 04 ob servation is in a s tate similar to the earlier 
02 reported in lGallo et al 1 (|2004l ). The high energy spectral 
drop has an energy of 7.311q" keV . The two models, par- 
tial covering and reflection, both give an eq ually good fit to 
the da ta above 3 keV . However, as noted in Revnold s et al.l 
(|2009l ) for MCG-6-3-15, an absorber that is responsible for 
the spectral drop is expected to produce fluorescent emis- 
sion at '-^6.4 keV . This is clearly not seen in the data, and 
might require a special geometry, where the absorber covers 
only a small part of the sky as seen by the source. The value 
of ~ 20 per cent covering fraction for the absorber at the 
source, is an upper limit, and although it is not as small 
as that reported in I Revnolds et al.l (|2009l ) , it points to the 
fact that a 'special geometry' explanation cannot hold for 
all AGN. Since, to reproduce the drop, most of the absorber 
must lie exactly in our line of sight with a covering fraction 
of 70 per cent, and a global covering fraction (i.e. as seen by 
the source) of about 20 per cent. 

The drop is more consistent with being the blue wing of 
a broad Ka emission line, originating very close to the black 
hole. The inner radius of the emitting region is ^ 1.3rg, and 
if we assume no radiation is emitted from within the radius 
of marginal stability, the black hole has to be an almost 
maximally-spinning Kerr black hole. A fit using kerrdisc 
model from Brenneman Revnoldsl (|2006l ) gives a spin pa- 
rameter a > 0.976 at 90 per cent confidence (where a is the 
dimensionless parameter a = cJ/ GM^ , with J b eing the an- 
gula r momentum of a black hole of mass M, see Miller et al.l 
I2OO9 for spin measurements using this method for stellar 
mass black holes). Fig. [24] shows the changes for the spin 
parameter. The parameter is well constrained to be > 0.97, 
and this is the case also if it is measured from the low flux 
observations (01 and t he first of 03). A lthough this might 
not be an exact value (iBrenneman R evnolds 2006, sub- 
mitted), what is clear from the shape of the line is that 
most of the radiation has to be emitted from within a few 
gravitational radii, and the black hole in 1H0707 appears to 




0.4 0.6 

Spin parameter a 

Figure 24. Ax^ versus spin parameter for a fit using kerrdisc. 
The black hole is rapidly spinning. The fit was made to 3.0-10.0 
keVband for the combined 04 observations. Inset: a zoom-in on 
the high spin values. The dashed line in both plots marks the 3cr 
limit. 



be rap idly spinning similar to what was found for MCG-6- 
30-15 ([Brennem an R evnoldsl [iooeh . 

This interpretation is further strengthened with the 
identification o f the feature at ~ 1 keV as the a broad iron L- 
sheh line ( Fabian et al.ll2009l l. As we have shown in Sec. 13.21 
(see Fig. [5]), the spectral complexities around 1 keV cannot 
be fitted with any absorption-only model, and there is a 
requirement for an emission component. The flux ratio be- 
tween the K and L features is 1:20, measured as the ratio of 
equivalent widths of two laor hues, is consistent with th e 
value expected from atomic physics (|Ross &: FabianI [20051 ). 

It is apparent from the models fitted to the high energy 
band that the spectrum requires high iron abundance. This 
is expected to imprint other signatures in the spectrum. In 
particular, if the apparent drop at 1 keV is interpreted as 
an iron L edge, then a deep absorption trough is expected to 
be produced by iron M-shell UTA. A reflection continuum on 
the other hand, accounts naturally for most of the features. 

The high iron abundance ('-^ 9x solar), if real, is not in- 
consistent with wh at is expected from line ratio modeling at 
other wavelengths ([Shemmer Netzerll20o3 : iHamann et al.l 
2002), and would imply that the host galaxy had gone 
through a phase of strong star formation. This seem to be 
consistent with studies showing that NLSl have st ronger 
star formation than normal Seyferts (Sani et 811120091 1. Also 
the Fell emission in also strong in these sources, and it ca n 
also be attributed to high abundance ([Collin &: Jo1vI[200q[ V 
The high abundance may also be due to a dense nuclear 
supernova activity. On the modeling side, non-solar values 
for elements other than iron can slightly affect the shape of 
the broad iron Ka and the abundances inferred from the fits 
(Rev nolds et al.lll995l ). 

Because the source is highly variable, fitting a time- 
averaged spectrum might not give the whole picture. Both 
flux and time-resolved spectroscopy indicate that the spec- 
trum is composed of two components, a highly variable com- 
ponent between 1-4 keV and a less variable component out- 
side the range. These two components are consistent with 
a PLC and RDC respectively, similar to the results from 
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other AGN (iFabian Vaughanl l2003l : iTavlor et all l2003l : 
IVaughan FabiaJ 20041 ) . In particular, the shape of the 
constant (or slowly varying) component dominating the low 
flux state matches very well the shape of a reflection spec- 
trum (see Fig. [T]). Variations in the power- law component 
drives most of of the variability (Fig. [15)), and this simple 
model reproduces the variability pattern very well (Fig. I12p . 

No absorption-only model provides a good fit to the 
spectrum, and these models can also be ruled out based on 
spectral variability. If the spectral drop at ^ 7 keV is due to 
an absorbing cloud in the line of sight, then the depth and 
shape of the drop imply the absorber only partially covers 
the source, which in turn means it has to be separate from 
any absorber causing the drop at ^ 1 keV . This however, 
is inconsistent with the variability pattern showing that the 
spectral features at ^ 7 and ^ 1 keV are varying in a similar 
way (Fig. [12]). 

The variability pattern, where a good correlation be- 
tween reflection and the direct continuum is seen at low 
fluxes and levels off at higher fluxes, is consistent with the 
predi c tions of the light b e nding model (|Fabian Vaughanl 
I2OO3I : iMiniutti k FabianI l2004h . Most of the luminosity 
emerges from within few gravitational radii, and strong grav- 
itational light bending must be taken into account. However, 
the model in its simplest form assumes that the intrinsic 
luminosity of the source is constant, which most likely is 
not the case. Also, assuming a single ionisation parameter 
throughout the disc might not be reaUstic too. These two 
factors may be responsible for producing the positive corre- 
lation in Fig. I16[ where an increase in the intrinsic luminos- 
ity of the primary component causes the disc ionisation to 
increase. The most likely scenario is that both light bending 
and intrinsic source variations are at work. 



6.1 Time Lag and Reverberation 

The lag spectrum is frequency dependent, and cannot be 
fitted with a single power-law. If however, the fit is only ap- 
plied to the central frequency bins, a good fit is found with 
t(/) oc /~°"^^. The slope cannot be constrained very well, 
and a power- law with an index of 1 is also consistent with 
the data. This is similar to (or at least consisten t with) what 
was found in NGC 7469 frapa dakis et al ]|200lh.MCG-6-30 - 
15 (|Vaughan et al.|[2Q 03). NGC_4051 (|McHardv et alll2004l ) 
and NGC 3783 ( Mar kowitzll 200 51. Several models have been 
suggested to explain the lag, the most suc cessful of which is 



the propagation o f accretio n fluctuations (lLYubarskiilll997l : 
iKotov et aLll200lL see also lArevalo Uttlevll2006l ). In this 



model, accretion rate variations are produced in different 
radii in the accretion disc, each with time-scale associated 
with the radius of origin. These fluctuations propagate in- 
ward and modulate the central X-ray emitting region. This 
seems to explain most of the timing properties including the 
PSD shape, its changes with energy and RMS-flux relation 
(|uttlevll2004l ). If additionally, energy-dependent emissivity 
profiles are assumed, with soft energies having flatter emis- 
sivity profiles, the dependency of the lag on frequency and 
energy can also explain the observations in galactic black 
holes and AGN. 

The lag is produced when large time-scale fluctuations 
propagating inward passes through the soft-emitting region 
first (larger radii), before reaching the hard-emitting region 



(|Kotov et al.l l200lh . This produces the positive lag below 
~5xlO~^Hz. The magnitude of the lag itself is dependent 
on the propagation speed and the distance between the soft- 
and hard-emitting regions. 

The slope of the lag spectrum appears to change 
at small temporal frequencies (large time-scales). If this 
turnover is real (i.e. the last point in Fig. [18]), it might be the 
effect of the size of the emitting region. The emissivity pro- 
files inferred from the spectral fitting to the iron line indicate 
that most of the radiation originates within few gravitational 
radii. This might pose an upper limit on the hard lags that 
can be observed. Interestingly, the maximum of the lag, at 
1 — 2 X 10~^ Hz, is about the same as the break frequency in 
the PSD, which was found to be 1.6lo e x 10"^ Hz, (when 
fitted with a broken power- law). The break in the PSD in 
this propagation model may be related to the time-scale of 
the fluctuations at the outer ra dius of the X-r ay emitting 
region ( Kotov et al. 2001; McHar dv et"all |2004| ) . The light 
crossing distance corresponding to ^ 150 s for a black hole 
mass of 5 X lO^M© is ~ 8 rg . The viscous time-scale at 
this radius in a standard disc (|Shakura Svunvaev|[l973h 
is tvisc 10^ s (assuming a viscosity parameter a = 0.1 and 
a scale height H/ R — 0.1)^ which is much higher than what 
is observed. If however, thermal time-scales are considered 
where tth = (i//i?)^tvisc (which might be the case if variabil- 
ity is caused by ionisation changes), then the lag peak and 
the PSD break would correspond to the thermal time-scale 
of the outer edge of an emitting region that is less than 8rg 
in size. 

At time- scales below about half an hour, the lag be- 
comes negative, where hard energies now lead softer en- 
ergies. This is the first si gnificant soft lag detected from 
black holes (|Fabian et al.lT2009.. V Other authors have re- 
ported soft lags befo re, but they wer e not significant ( e.g. 
iMcHardv et"al1l2007 l for Ark 564 and lGallo et"alll2004l for 
IRAS 13224-3809). If the positive lag is interpreted as being 
due to Comptonisation, where the variability is seen first in 
the soft band (disc component in the form of reflection) be- 
fore the hard band (Comptonised coronal emission), then 
negative soft lag above ~ 6 x 10~^ Hz is simply in the 
wrong direction. More sophistica ted models may be tuned to 
produce soft lags. In particular, iMalzac k, JourdainI (|2000h 
studied the temporal evolution of flares in a disc-corona sys- 
tem. They found that, if a flare originates in the corona (e.g. 
by magnetic reconnection) with the perturbation time-scale 
of the order of few corona light crossing time, then a soft lag 
is expected, and is caused by a feedback mechanism that 
cools the corona down. Fluctuations are seen in harder en- 
ergies, before reprocessed soft photons from the disc cool the 
corona and softens the spectrum. In this picture, the whole 
spectrum originates in the corona. This however, is incon- 
sistent with the spectral variability patterns seen in Sec. H] 
There is an abrupt change in variability properties between 
the bands above and below ^ 1 keV (see also RMS spec- 
trum in Fig. 8 of the supplementary material in Fab ian et all 
12009), which argues for two separate components above and 
below 1 keV. 

In a reflection scenario, the soft lag is naturally ex- 
pected. The variability is seen first in the power- law dom- 
inated component, emitted from the corona. Some of the 
radiation is reflected off the optically thick disc. This pro- 
duces a lag comparable to the light travel time between the 
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the corona and the reflector. If we assume most of the ra- 
diation is emitted at ^ 2rg, imphed from the inner radius 
and emissivity index in the spectral fitting, then interpret- 
ing the 30 s lag as the light crossing time yields a black hole 
mass of M ~ 2 X lO^M©, which is consistent with the un- 
certain mass of this black hole quoted in the literature (e.g. 
IZhou Waii3l2005l l. 

Ark 564 shows a similar la^ spectrum to the one pre- 
sented here (McHard v et al.ll2Q0/t ). with a difference in the 
magnitude of the lag. The ratio of hard to soft lags in Ark 
564 is - 700/10 = 70, compared to - 5 in 1H0707. This 
appears to be related to the compactness of the emission re- 
gion and the strength of the reflection. 1H0707 has a strong 
reflection component associated with a compact emission 
region with strong light-bending effects. However Ark 564 
has weaker reflection which suggests a more extended emit- 
ting region. Since thermal and viscous time-scales scale as 
(r/r^)*^/^, the hard lags due to changes in the accretion flow 
in this region are larger compared to the light-travel time. 

Unlike other well-studied AGN (e.g. MCG-6-30-15 and 
Mrk 766), the spectrum of 1H0707 does not show any signs 
of warm absorption, giving an uninterrupted view to the in- 
ner regions of the accretion system. 1H0707 h as one of th e 
strongest soft excesses among other NLSl fe.g. lGall"3l2006h . 
That, in our interpretation, is simply because of the high 
iron abundance making the iron L complex st rong, which 
might also be linked to the conclusion of Gallo (2006), that 
most NLSl with complex spectral features show high val- 
ues of Fe ii/H/3. Other sources that are similar to 1H0707 
include IRAS 13224-3809, which shows clearly the iron L 
broad feature, however the current XMM observation is not 
long enough to check if there is any soft lags (Ponti et al. 
2009, submitted). 

In summary, although different models may be made to 
explain some of the properties we see independently, most 
of them fail to explain all of what we see. In this work, we 
have explored the data using different methods, and shown 
that a model in which the contribution of X-ray reflection 
off optically thick material is significant, provides a natural 
explanation to the spectrum, spectral variability, as well as 
the soft lag seen for the first time in 1H0707. 
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